BCS— BEC crossover in dense relativistic matter: Collective excitations 
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We study the relativistic BCS-BEC crossover within a class of Nambu-Jona-Lasinio type models, 
including arbitrary Lorentz-scalar pairing channels. Using the mean-field approximation we investi- 
gate spectral properties of the collective bosonic excitations in the superfluid phase, with particular 
attention to the Nambu-Goldstone bosons of the broken symmetry. This is a first step towards a 
systematic improvement of the mean-field approximation by including the fluctuation effects. The 
general results are illustrated on pairing in dense two-flavor quark matter — the two-flavor color 
superconductor. 
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I. INTRODUCTION 

It was noted long ago that Bardeen-Cooper-Schrieffer 
(BCS) superconductivity of correlated Cooper pairs and 
Bose-Einstein condensation (BEC) of tightly bound 
difermion molecules are merely two sides of the same 
coin. Eagles EI and Leggett 0] were perhaps the first 
to observe that a many-fermion system can display both 
the BCS and the BEC behavior, and these two regimes 
are connected by a smooth crossover as the strength of 
the attractive interaction between the fermions is varied. 
The crossover can be described at least qualitatively well 
by the conventional mean-field approximation provided 
one determines the chemical potential self-consistenly by 
fixing the overall particle density. 

However, it was only recently that the BCS-BEC 
crossover could be realized experimentally using ultra- 
cold atomic Fermi gases [!, 0, H, II , where the strength 
of the interaction can be tuned by external magnetic 
field and thereby, using the Feshbach resonance, cover 
the whole range between the BCS and BEC limits. The 
dilute atomic gases represent particularly clean and well 
controlled examples of interacting many-body systems. 
It is believed that they may provide an ideal theoretical 
playground for development of many-body techniques, 
that would be later applied to more complicated strongly 
correlated systems. 

The mean-field theory of Eagles and Leggett was for 
the first time improved by Nozieres and Schmitt-Rink 
Q who included the contribution of the pair fluctua- 
tions to the total particle density. In the past decade, 
various extensions of the mean-field approximation have 
been proposed, including boson fcrmion models 
well as purely fermionic self-consistent schemes [9, 
An extensive list of references may be found in the re- 
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views [ll|, |l2| . The most recently developed techniques 
are compared in El . Simultaneously to the analytic ap- 
proximation schemes, the interacting Fermi gas is also 
being studied using numerical Monte Carlo simulations 

El si. 

In high-energy physics, the analogy with BCS super- 
conductivity was used long ago by Nambu and Jona- 
Lasinio [la , Il7j to propose a model for dynamical chi- 
ral symmetry breaking. Later, Cooper pairing of two 
fermions near their Fermi surface found an application 
in the physics of dense nuclear/quark matter, giving rise 
in particular to the phenomenon of color superconduc- 
tivity El El; 

see |20] for the most recent review. The 
crossover physics is relevant for these high-energy sys- 
tems because of the strong-coupling nature of quantum 
chromodynamics. Concretely, the so-called pseudogap 
phase as a precursor to superfluidity of tightly bound 
Cooper pairs was investigated in models of chiral sym- 
metry breaking [2l|, E2, EH HI as we ^ as color super- 
conductivity (25| • The structural change of Cooper pairs 
with increasing coupling strength and the possibility of 
their BEC were studied in [H, [27|, El, El ■ The spectrum 
of diquarks in various color-superconducting phases was 
investigated by Ebert et al. [30, [3l| ■ 

The actual crossover between the BCS and BEC 
regimes in dense quark matter in dependence on the cou- 
plin g streng th has started to be studied only recently 
[33 . l3a . [3J|- Perhaps the first attempts to go beyond 
the mean-field approximation have been made in Refs. 
HEEEE3- In particular, Abuki (H has conducted an 
extensive study of fluctuation effects in the normal phase, 
using the Gaussian approximation. The goal of this pa- 
per is to investigate the spectrum of the bosonic collective 
modes below the critical temperature. Thereby we con- 
struct a framework for description of fluctuation effects 
in the superfluid phase. 

The plan of the paper is following. In the next sec- 
tion, we derive the general formulas valid for a class 
of models of the Nambu-Jona-Lasinio (NJL) type, that 
describe (even-parity) spin-zero pairin g in relativistic 
fermion matter. The model of Abuki [391, which is a 
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special case, is then investigated in detail in the follow- 
ing section. In particular, we study the spectral prop- 
erties of the bosonic collective modes in the mean-field 
approximation. Possible extensions of the present work 
are discussed in the conclusions. 



II. NJL MODELS OF PAIRING 

A. Model definition 

We consider the general class of models defined by the 
Lagrangian 



C=r(,il) f ,,,„- ,»lr + ^5^(^ T r ^)(^f o ^ r ), (1) 



where T a = jqT^jq. The spinor ip denotes a set of 
fermion flavors with common mass m and chemical po- 
tential /i, while T a are matrices that furnish an an- 
tisymmetric tensor representation under the symmetry 
transformations acting on ip. In general, they act on 
both Dirac and flavor indices of ip. Even though the 
Lagrangian |T]) covers, in a simplified manner, a num- 
ber of systems investigated in literature [38|], at this 
stage we completely neglect correlations in the fermion- 
antifermion channel that may lead to dynamical breaking 
of the chiral symmetry. 

This model is analyzed using the standard Hubbard- 
Stratonovich transformation. We start with the defini- 
tion of the Nambu spinor 

charge conjugate field, ip c = Cip T , the interaction term 
may be rewritten as § |^> C P ^>| , where P a = C~ x T a . 
It is removed from the Lagrangian by introducing new 
auxiliary scalar fields A a and adding 
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The new, semibosonized Lagrangian reads 
r . — ——A* A + ivfrD -1 ^ 

•'-semi — q L - x a L - i a T 2 



where 



p_! _ (0 + 7oM - m A a P Q 



A*P Q i0 - 7o// -ml ' 



(2) 



is the fermion propagator in the coordinate space. Upon 
integrating out the fermions, the thermodynamic poten- 
tial £1 is given by 

e-" n = ydA a dA;exp(-S off [A,A*]), 

/ dr /" rf 3 x|A a (x,r)| 2 -ilogdet2?- 1 . 

(3) 



This functional integral, of course, cannot be evaluated 
exactly. Note that the bosonized action S e s includes the 
effect of fermion loops to all orders, due to the presence of 
the determinant of the inverse fermion propagator T)^ 1 . 
In other words, it contains all Fcynman graphs with the 
scalar fields as external lines. It is thus to be interpreted 
as a classical action for the scalars A a , A*, and the above 
functional integral as the generating functional of their 
Green's functions (at zero source). This can be deter- 
mined from the quantum effective action associated to 
the classical action S e g, being the sum of all one-boson- 
irreducible graphs. 

Since we assume that the leading nonpcrturbative be- 
havior, which gives rise to the pairing, is borne by the 
fermion loops that have already been resummed to all 
orders, the expansion in the number of bosonic loops pro- 
vides a natural approximation scheme for the calculation 
of the thermodynamic potential. The tree level corre- 
sponds to simply setting /3£Imf — S e s'. This is the stan- 
dard mean-field approximation. At one loop, the ther- 
modynamic potential would be given by 



= /3^mf + - log det S, 



(4) 



where 



E ab {X-Y) = 



S 2 S cff 



( xff * 

A* A* 



6$l(X)6$ b (Y) \Xab 



AA ■ 

^A ) , (5) 

Xab 



is the inverse propagator of the bosonic modes; $ is the 
scalar doublet field, <I> Q = (A a , A*) T . Note that in the 
normal phase, Eq. ((4]) reduces to the Gaussian approx- 
imation adopted by Abuki (3r| . Carrying out the indi- 
cated differentiation and setting the fields A a equal to 
their vacuum expectation values (which we hereafter re- 
fer to by the same symbol), we arrive at the expressions 
for the normal and anomalous parts of the inverse prop- 
agators, 



n ^ ' 

x Tr [p a P**(iw n ,k)nP^(»w B - ifijv.k- p)j , (6) 

x IV [PJ}**(ko„,k)P h V**(iw„ - m„,k - p)] , (7) 

where ip and ip in the superscripts to T> denote the re- 
spective matrix elements of the matrix propagator |2|). 
The remaining two entries of the matrix ([5]) follow from 
the general relations, 

x£ a (",p) = x& a *(-",-p), 

AA/ \ AA/ \ 
Xab (W,P) =Xba (-<«>, -P), 

valid for an arbitrary complex frequency tu. 
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B. Fermion propagator 

In order to be actually able to invert the matrix "D -1 
@ , we now make a specific assumption about the fermion 
pairing pattern. We assume that the Cooper pairs are 
Lorentz scalars, i.e., the matrices T a have the form T a = 
C"/c,Q a , where Q a act just on the internal symmetry in- 
dices. The order parameters A a now enter the fermion 
propagator in terms of the expression M = A a Q a . Note 
that Pauli principle requires T a to be antisymmetric so 
that both Q a and M have to be symmetric matrices. 

With the above assumption, the denominators of the 
matrix fermion propagator depend only on the combi- 
nations MM^ and M'M. Being Hermitian and positive 
definite, the matrix MM' may be spectrally decomposed 
as [H 

MM* = J2A 2 r V r , 

r 

where V T is a set of projectors on the respective eigen- 
vectors, and Af, are the real positive eigenvalues whose 
relation to A a will be clarified later. Using the symmetry 
of M, we find the equivalent decomposition, 

m^m = ^A 2 r v;. 



The parameters A r play the role of gaps in the 
(fermionic) quasiparticle dispersion relations, which are 
found as poles in the fermion propagator D, 



with = e k + efi, e k = vl? 



(8) 



Explicit expressions for the matrix elements of the 



fermion propagator, obtained by inverting Eq. p|). read, 

ioj n - e£ k 



r e— ± 



r e— ± 



(9) 



r e— ± 



V r M, 



A k 7s 



r e— ± 



The definition and basic properties of the standard en- 
ergy projectors Ajj. are given in Appendix I A H 



C. Propagator of the collective modes 

Eqs. ©,0, and @ show that the calculation of the 
inverse propagator, H, of the bosonic collective excita- 
tions can be split into three independent steps. First, 
the trace over Dirac indices is always of the form 



Tr D (A£A/) = l + e/- 



kq 



£k£c 



Second step is the trace over the flavor indices, TrF. This 
will be discussed shortly. Last, the summation over the 
fermionic Matsubara frequencies uj n gives rise to func- 
tions i^/(k, q;iQ;v) and J®/(k, q; if2jv) defined in Ap- 
pendix E2J 

Putting all the pieces together, the normal and anoma- 
lous correlation functions, x AA and x AA , are given by 
the formulas 



xT'(^n, P ) = ^Sa b + lj2TT F (Q a V r QtP:)J2j ^Tr D (A^A(_ p )4 e /(k,k-p;z^), (10) 
X^ A (*n W ,p) = ^Tr F (Q a V r MQ b V s M)Y^ / 7^ Tr D (A k A£_ p ) ^/(k, k - p; iSl N ). (11) 

r,s e,f ^ ' 

In particular, in the normal phase all A's are equal to zero so that x AA = 0. Moreover, we have 
s T I F(Qa'PrQl'Ps ) = Trp(QaQb) = N5 a b, which can always be enforced by a suitable choice of basis of the 
flavor algebra. N is a normalization factor which, in principle, can differ for different pairing channels (or, different 
irreducible representations of the flavor symmetry). In the following, we will assume that the fermions pair in a single 
channel, but the generalization is obvious and straightforward. After this remark, we can write down the correlation 
function in the normal phase as a special case of Eq. (fTU|) , 

1 - 2/(e k - fi) 

AT + 2fl — £ k — Ck-p. 

/(ek + ^)-/(ek- q -^) ] 
ifl N + 2/ji + e k - e k _ p J ' 



N 

x(ift N ,p) = — + — 



d 3 k 

(2tt) 3 



k(k-p) 



CkEk-r 



2/(e k + /x)-l 



+2 



2fi 



' £k + £k-p 

I 2 + k ■ (k - 

£k£k-p 
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which is equivalent to Eq. (15) in Ref. [361 ]. 

Full information about the spectrum of bosonic collec- 
tive states may be obtained by investigating the spec- 
tral density, p(u>,p), associated to the matrix correlation 
function S, defined by 



J — oo " 



du) p ab (ui,p) 



(12) 



In general, we would like to use the relation 

p ab (u,p) =lm[E,- 1 (uj + iS,p)] ab: (13) 

in order to extract the spectral information from the 
imaginary part of the pair correlation function, analyti- 
cally continued to the real axis. However, this should be 
done with some care for Eq. (fT3| holds only when the 
spectral density p a b(u>,p) defined by Eq. (|12[) is real. 

D. Mean- field analysis 

Assuming that in thermodynamic equilibrium the sys- 
tem relaxes to a homogeneous state, the mean-field ther- 
modynamic potential ^mf is given by Eq. ([3]), 



MF 



V 



G 



(14) 



where we subtracted the energy of the vacuum in order 
that f^MF expresses the sheer effect of finite temperature 
and density, and the Cooper pairing. 

It should be noted that Eq. (fT4")) embodies two different 
sets of gap parameters. First, the A a 's, which carry the 
label of the Cooper pair representation of the symmetry, 
and represent the order parameters for symmetry break- 
ing. Second, the A r 's, which determine the gaps in the 
fermion excitation spectrum. While the A a 's transform 
as an antisymmetric tensor under the symmetry trans- 
formations acting on tp and may in general be complex, 
the A r 's are by construction invariants of the symmetry 
and are always real. 

These two sets of parameters may be related with the 
help of the expression A a A*Qt Q b = MM^ = £ r A 2 -p r . 
By taking the trace we find 



iV^|A Q | 2 = ^A^. 



This clarifies the physical content of the normalization 
factor N: It counts the number of fermion degrees of 
freedom participating in the formation of a particular 
Cooper pair. Plugging this result back into Eq. (fl"4|) . 
we immediately arrive at the standard mean-field gap 



equation and the equation for the particle number density 
n, 



e— ± 



d 3 k l-2/(££ r ) 



2E\ 



kr 



r e— ± 



(2tt) 3 Et r 



(15) 
(16) 



The total number density n is fixed by the value of the 
Fermi momentum fep, as n = fcp/ (37r 2 ) times the number 
of fermion flavors. 



E. Renormalization 

The momentum integral in Eq. (|15p (as well as in other 
equations depending explicitly on the coupling) is badly 
divergent. This can be taken care of by trading the bare 
coupling G for the physical s-wave scattering length at 
zero temperature and density, or the renormalized cou- 
pling Gr, defined by [36[ 



1 

G 



1 

Gr 



N 

y 



d 3 k 

(2tt) 3 



As opposed to the case of nonrelativistic fermion mat- 
ter, the remaining momentum integration is still diver- 
gent, but only mildly, logarithmically. In numerical com- 
putations, it is regulated with a sharp three-momentum 
cutoff A. All equations of the theory may be rewritten 
in such a way that the cancelation of the leading diver- 
gences is manifest. The resulting formulas are useful for 
the numerical implementation, but otherwise are rather 
cumbersome. Just for illustration, we show here the gap 
equation and the equation for number density, 



A r 



1 



A r NG 



*E 

e=± 



d 3 k 

(2ir) 3 



1 2 

k e k 



E 



kr 



2p 2 



^kr-^kr i^k 



EE 

r e— ± 



d 3 k 
/"A 2 . 



-2e-^f(E£ r ) 



Et E7 e (Ei +EZ e ) 

kr kr \ kr 1 kr / 



E\ 



kr 



E, 



kr 



£k 



The procedure outlined above is sufficient to renor- 
malize the mean-field thermodynamic potential and its 
derivates, the gap and number equations. On the other 
hand, one must be more careful when calculating the 
correlator S a t> ([3]), or in general any fermion loop with 
nonzero external momentum. The point is that the mo- 
mentum assignment for the internal fermion lines is not 
completely fixed by momentum conservation in the in- 
teraction vertices; it is unique up to a shift of the inte- 
gration variable. Such a shift, however, is not permitted 
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in an integral with stronger than logarithmic divergence 
(50| . Therefore, we are led to consider the momentum 
assignment to the fermion propagators in Eqs. ([6]) and 
([7]) as a part of the definition of the model. As already 
explained above, the integrals are then evaluated with a 
sharp cutoff on the three-momentum k. 

Let us finally remark that one might at first sight think 
it would be natural to impose the cutoff directly on the 
arguments of the fermion propagators in the loop, i.e., 
on both k and k — p in Eqs. ^ and (J7J); this would 
roughly correspond to discretizing the coordinate space. 
Nevertheless, this would destroy the expected low-energy 
dynamics of our model, based on the Nambu-Goldstone 
(NG) modes of the spontaneously broken symmetry - 
the inverse propagator ((6J would develop a linearly di- 
vergent piece, proportional to |p|. 

III. TWO-FLAVOR COLOR 
SUPERCONDUCTOR 

Now we apply the general formulas derived in the pre- 
vious section to the particular model studied by Abuki 
[36| . It describes color superconductivity in quark mat- 
ter consisting of two quark flavors, which have for sim- 
plicity equal masses and chemical potentials. The field 
ip now carries the index of the fundamental representa- 
tion of the symmetry group SU(3) x SU(2). Since the 



weak-coupling studies of quark matter indicate attrac- 
tion between two quarks in the color-antitriplet chan- 
nel, we assume that the matrices Q a have the structure 
(Qa)b c — £ lJ £abc, where i,j from now on denote the fla- 
vor SU(2) indices and a,b,c the color SU(3) ones. The 
order parameter A Q thus transforms as an SU(2)-singlct 
and an SU(3)-antitriplet. As usual it is chosen to be 
real and to point in the third (anti-blue) direction in the 
color space. In the following, we write simply A instead 
of A 3 ; this will help us distinguish the constant order 
parameter from the fluctuation fields without having to 
introduce further notation. 

For this particular symmetry structure and the choice 
of vacuum, we find N = 4 and = A 2 diag(1, 1,0) 

in the color space (in the flavor space, it is simply the 
unit matrix). So in this case, the projectors V r sim- 
ply project on the individual colors. The red and green 
quarks, participating in the pairing, are gapped, while 
the blue quarks remain gapless. As is well known, the 
global SU(3) x U(l) symmetry of the bosonized action 
([3]) is broken by the Cooper pair condensate A down to 
SU(2) x U(I)q, operating exclusively on A x and A 2 [5l| . 
The unbroken symmetry implies degeneracy between Ai 
and A2 . The pair correlation matrix S a b is then diagonal 
and Sn = S 2 2- Moreover, En has only the normal part, 
thanks to the existence of the unbroken U(1)q charge. 
Evaluating explicitly the color-flavor traces in Eqs. (fTU)) 
and pip , we obtain the expressions 



*(iO iV ,p) = i + 2^ J 

e,/ 

X^ 3 A (^,P) = -2A 2 ^ J 



d 3 k 

(2tt) 3 

d 3 k 



1 + e/ 



m 2 + k • (k — p) 



(2tt) 3 

d 3 k 

(2tt) 3 



1 + e/ 



EkEk-p 

to 2 + k • (k — p) 



1 + e/ 



CkEk-p 

in 2 + k • (k - p) 



£kek-p 



[7 2 t(k, k - p; in N ) + 7 3 e |(k, k - p; iil N )], 



A £ 2( k ; k - p;*^jv), 
Jl£(k,k-p;in N ). 



Note that the quark loop contributing to S33 is symmet- 
ric: Both propagators correspond to gapped quarks. On 
the other hand, the loop in En consists of one gapped 
and one ungapped propagator and the general formula 
(I10p requires it to be symmetrized. 



A. Nambu— Goldstone bosons 

Spontaneous symmetry breaking leads to the existence 
of NG bosons, soft fluctuations of the order parameter. 
Vacuum expectation value of a color antitriplet breaks 
five generators of the color SU(3) so that we should ex- 
pect five NG bosons. However, in our model where the 
color symmetry is just global unlike in quantum chromo- 



dynamics, the Noether charge corresponding to the Gell- 
Mann matrix Ag acquires nonzero density. This implies 
that there are only three NG bosons, two of which have 
quadratic dispersion law at low momentum, being the 
so-called type-II NG bosons [H Ell, 53] . These couple to 
the two pairs of generators, (A^As) and (Ag,Ay), whose 
commutator contains a As piece and hence has nonzero 
density. Therefore, they will be found in the spectrum 
of excitations of Ai and A 2 . The remaining, type-I NG 
boson couples to the spontaneously broken linear combi- 
nation of A3 and As and will be found as an excitation of 
A 3 . 

At zero temperature, the presence of NG modes in our 
model may easily be demonstrated analytically. We set 
p = 0, use the gap equation (TH)]) , and analytically con- 



tinue the pair correlation functions to real frequency (or, 
more precisely, to the vicinity of the real axis), 



1 



1 



-ST) 



d 3 k 1 uj 2 - 2A 2 - 2ew^ 



e=± 



Xssf ( W ' )- 2 E / ( 27r )3 £e ^2 _ (2^)2 



d 3 k 1 1 

(2^^-(2££)' 



[We have dropped the index r so that is now defined 
by Eq. ([5]) with the gap A inserted.] We can immediately 
see that Xii A (^j 0) ~ > as ui — > 0, which proves the 
existence of NG excitations coupled to Ai and A2. These 
are therefore interpreted as massless bound states of one 
gapped and one gapless quark (green and blue, and red 
and blue, respectively). 

The correlation functions of A3 do not go to zero at 
vanishing frequency However, here the propagator is 
given by inverting the matrix S ([5]). Since in the limit 
u> — > we have 

x£f(o,o) = -x& A (o,o), 

the determinant of S will be zero and the existence of a 
NG mode is thus proved. 

At zero temperature, we can of course go beyond the 
mere demonstration of the existence of NG bosons, and 
determine their low-momentum dispersion relation. This 
may be done in the standard manner, by expanding the 
correlation functions in powers of momentum. Let us 
just remark that at nonzero temperature but still in the 
superfluid phase, there are kinematical regions where the 
sharp NG peak in the pair correlation spectrum is ob- 
scured by Landau damping — the emission or absorp- 
tion of the NG boson by a quark. In such regions the 
momentum expansion of the correlation functions is not 
possible [HI]. 

Here we concentrate on the doublet of type-II NG 
bosons coupled to Ai : 2, aiming to investigate the effect of 
the expected quadratic onset of their dispersion relation. 
At low frequency and momentum, the inverse propagator 

All 



AA is expanded as 



Xu A >,P) 



bp 2 



(17) 



The coefficients a, b determine the low-momentum dis- 
persion relation of the NG boson by w(p) = p 2 b/a. They 
are given by the explicit expressions 



a = A^E 



e=± 



rf 3 k 



H-sgn(e^) 




G /G_ 

c R 

FIG. 1: The critical temperature (solid) and chemical poten- 
tial (dashed) in units of the Fermi energy, as a function of the 
inverse coupling. Three different regimes are distinguished: 



BCS (Gc/Gr 
(G c /G r > 1). 



< 0), BEC (0 < Gc/Gr < 1), and RBEC 



b = 0{n — to) 



(p 2 - TO 2 ) 3 / 2 

3tt 2 AV 



1 - 



34 



1 



e=± 

s g n ?k 



(2tt) 3 



— e J 



2(^) 2 



2k 2 1 



k^ 
2 A 2 

e\5 



Note that a is actually equal to 1/ A 2 times the difference 
of the density of a gapped and a gapless color [cf. Eq. 
(fl"6|) . the factor 2 arises from the existence of two quark 
flavors for each color] . This in turn is proportional to the 
density of Ag as it should, for it is this non-Abelian charge 
which, according to general theorems [12, E3] gives rise 
to type-II NG bosons. 



B. Numerical results 

In this section, we present numerical results that illus- 
trate the general conclusions made above. In order to 
achieve concrete numbers, we use the same set of param- 
eters as Abuki [36| so that our results in the superfluid 
phase can be matched to his for the normal phase. In 
particular, we choose the Fermi momentum = 0.2m 
and a cutoff to regulate the logarithmic-divergent inte- 
grals as A = 5m. The renormalized coupling Gr is rep- 
resented in units of the critical coupling G c , at which the 
mass of the diquark molecule in vacuum goes to zero. 
This signifies the instability of the vacuum itself with 
respect to pair formation, and the onset of relativistic 
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FIG. 2: Zero-temperature density of quarks of a single paired 
(solid) and unpaired (dashed) color at zero temperature, ex- 
pressed as a fraction of the fixed total density. Recall that 
there are two paired colors (red and green) and one unpaired 
one (blue). 



Bose-Einstein condensation (RBEC). 



1. Phase diagram 

The dependence of the critical temperature on the 
renormalized coupling is displayed in Fig. [1] It is in- 
structive to compare this mean-fie ld p lot to a similar one 
achieved by Abuki and Nishida [35l |36| . who included 
the effects of Gaussian fluctuations. (To be precise, they 
combine the mean-field Thouless criterion [45| with the 
equation for particle number that includes the contribu- 
tion of bosonic quasiparticles.) It is clearly seen that as 
expected, the mean-field approximation is reliable in the 
BCS regime (negative coupling), while it fails by about 
an order of magnitude on the BEC side. The reason is, 
of course, that the mean-field thermodynamic potential 
includes only the contribution of the fcrmions, whereas in 
the BEC regime the total number density is dominated by 
the bosons. In the RBEC limit (that is for G c /Gr ^ 1), 
both bosons and fermions are excited significantly, and 
the reliability of the mean-field approximation again im- 
proves. 

In Fig. [2] we show the contributions of the paired and 
unpaired quarks to the total density at zero temperature. 
In the BCS limit the unpaired blue quarks form a sharp 
Fermi sea, while the paired red and green quarks occupy a 
Fermi sea, smeared about the Fermi surface at the scale 
of the gap. However, at G c /Gr rs 0.08 the chemical 
potential decreases below the fermion mass and the Fermi 
sea can no longer exist. This marks the onset of the BEC 
regime. From this point on, the whole density at zero 
temperature is provided by the condensate of the bound 
diquark molecules, made of red and green quarks. The 
blue quarks can only be excited thermally. 




_l Ll_l i I i i i I i i i I i i i L 
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FIG. 3: Normal component of the spectral density matrix 
at zero temperature and momentum for G c /Gr = 0.5. The 
large-scale plot applies equally well to pn and P33. In the 
inset we show the tiny contribution of the particle-particle 
continuum at negative frequencies to P33, which is absent in 
pu- 



2. Bosonic spectrum at zero temperature and momentum 

As already remarked before, the correlation function 
S a b is diagonal in the color-flavor space, owing to the 
choice of the order parameter A to point in the third di- 
rection in the anticolor space. Moreover, since we choose 
A real, the spectral density p a b defined by Eq. ([12]) is 
real and we can use Eq. (fT5|) to calculate it. 

In the normal phase, the matrix S is simply propor- 
tional to the unit matrix and there is a single spectral 
density, which has been calculated at the critical tem- 
perature by Abuki 36]. In the superfluid phase, we have 
one spectral density for a = b = 1,2 and a symmetric 
2x2 matrix for a = b = 3. In Fig. [3] we give a sample 
numerical calculation of the normal parts of these spec- 
tral densities for G c /Gr = 0.5. 

Several comments are in order here. First, in the nor- 
mal phase an isolated bound state (corresponding to a 
diquark or antidiquark molecule) exists only on the BEC 
side of the crossover. In the superfluid phase, the NG 
bound state is present for all values of the coupling, as 
guaranteed by the Goldstone theorem. On the other 
hand, the antiboson pole appears only on the BEC side 
of the crossover, just as in the normal phase. It should 
be stressed that the positions of the antiboson poles in 
pn and P33 differ, simply because one of the bosons is 
composed of a gapped and an ungapped (anti)quark, 
while the other one of two gapped (anti)quarks. In- 
deed, the pole in p 33 occurs at a smaller (negative) u>; 
the corresponding boson is heavier. Nevertheless, for 
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FIG. 4: Absolute value of the position of the antiboson pole 
in (On (lower solid line) and P33 (upper solid line) at zero 
temperature and momentum. For comparison, we also show 
the values of 2 A (dashed) and 4p (dotted). 

Gc/Gr = 0.5 the difference is so small that it cannot 
be seen in Fig. [3J 

Apart from the presence of the isolated antiboson pole 
on the BEC side of the crossover (and its absence on 
the BCS side), the qualitative appearance of the spec- 
trum does not change as the coupling is varied. As noted 
above, the antiboson pole has different positions in p\\ 
and P33. The same, of course, holds for the onset of the 
continuum in these spectra. As we proceed from the BCS 
to the BEC to the RBEC regime, the gap in the spec- 
trum of red and green quarks increases, and the difference 
between pn and ,033 becomes more pronounced. 

In Fig. 0] we show the coupling-dependence of the 
two antiboson masses. In the BEC regime, they are 
almost equal and essentially determined by the chemi- 
cal potential. However, this changes dramatically in the 
RBEC regime. While the pole in P33 appears near the 
two-antiparticle continuum, now dominated by the large 
value of the gap, the antiboson pole in p\\ is still driven 
by the chemical potential. We checked numerically that 
as the temperature is increased, the mass of the anti- 
boson in P33 follows (twice) the value of the gap until 
very close to the phase transition, where it merges with 
the mass of the antiboson in p\\. This is in agreement 
with Ref. [36[ where it was shown that in the normal 
phase, both masses (being equal) are determined by the 
chemical potential. 

Another fact to note is the order-of-magnitude sup- 
pression of the anomalous spectral density P33 . This 
is due to the fact that, with the exception of the RBEC 
limit, the gap is much smaller than other scales in the the- 
ory such as the Fermi energy or chemical potential. This 



tiny off-diagonal (in the Nambu space) element of the cor- 
relation function matrix 5 then induces mixing between 
the bound boson and antiboson states, or in other words, 
breaks the conservation of the baryon number carried 
by the quarks. Consequently, the spectral density /o|^ A 
displays, besides the expected two-particle continuum at 
positive frequency, also an analogous continuum at neg- 
ative frequency, yet with much smaller spectral weight. 
In pf[ A , this continuum is absent due to the existence 
of the conserved charge of the unbroken U(1)q. 

3. Bosonic spectrum at nonzero temperature and 
momentum 

At nonzero temperature and momentum the spec- 
trum becomes more rich. The bound-state peaks acquire 
nonzero (even though tiny) width and the Landau damp- 
ing appears. Having in mind that the low-energy dynam- 
ics of the system in the superfluid phase is dominated 
by the NG bosons, we display in Fig. [5] the positive- 
frequency parts of the spectral density pn; each panel 
collects the spectra for several values of the momentum. 
The NG excitations are manifested by isolated peaks be- 
low the two-particle continuum. Comparison of spectra 
for temperatures close to zero and to the phase transi- 
tion reveals the temperature effects: The overall increase 
of the background due to thermal excitations, and the 
Landau damping, especially strong in the RBEC regime. 
We can also see how with increasing momentum the NG 
boson peaks get closer to the continuum; at a certain 
point they will reach the continuum and disappear from 
the spectrum as isolated poles. 



4- NG boson dispersion relations 

In addition to the full spectra, we also studied the dis- 
persion relations of the NG bosons at zero temperature. 
These can be extracted from the spectral densities and, 
in case of the type-II NG bosons, compared to the an- 
alytical low- momentum prediction from Eq. (|17p . The 
result is shown in Fig. [51 

General consequences of spontaneous s ym metry break- 
ing in presence of nonzero charge density [42| require that 
the dispersion relation of the NG boson in pn and P33 
is, at low momentum, quadratic and linear, respectively. 
The scale that defines this low-momentum region is typ- 
ically set by the symmetry breaking. Above this scale, 
the NG boson dispersion relation is no longer fixed by 
the general properties of symmetry breaking, but rather 
is determined by the detailed dynamics of the system un- 
der consideration. In some cases, the NG boson may even 
cease to exist at high enough momentum; in our model, 
this is reflected by its disappearance in the continuum. 

We display the dispersion relations in a logarithmic 
scale in order to exhibit their power-law behavior. Ap- 
parently, the low-momentum region governed by the 




FIG. 5: Positive-frequency part of the pn spectral density at nonzero temperature and momentum. First column: BCS regime, 
G c /Gr = —0.35; second column: BEC regime, G c /Gr = 0.5; third column: RBEC regime, G c /Gr = 1.35. The upper 
line shows spectra at low temperature, T = 0.05T C , the lower line near the phase transition, T — 0.95T C . The four curves 
in each spectrum correspond respectively to |p| = 0.5m (solid), |p| = 1.0m (dashed), |p| = 1.5m (dotted), and |p| = 2.0m 
(dash-dotted). 




p/m p/m p/m 



FIG. 6: Dispersion relations of the NG bosons. First plot: BCS regime, G c /Gr — — 0.35; second plot: BEC regime, G c /Gr — 
0.5; third plot: RBEC regime, G c /Gr = 1.35. In all graphs, the dispersion relation of the type-II NG boson in pn (solid), 
its low- momentum prediction (dotted), as well as the dispersion relation of the type-I NG boson in p33 (dashed) are shown. 
The hatched-gray and solid-gray regions denote the continuum in pn and P33, respectively. The dashed gray line indicates the 
value of the gap (BCS, BEC), or the chemical potential (RBEC). 



broken symmetry is indeed defined by the symmetry- 
breaking scale, i.e., the gap. The only exception is the 
RBEC limit where the predicted quadratic dispersion of 
the type-II NG boson is observed only below the scale of 
the chemical potential, which is much smaller than the 
gap. This is because above the chemical potential the 
particles and antiparticles become nearly degenerate and 



all effects of finite density are strongly suppressed. 

In the BEC regime, we can even distinguish three re- 
gions of momentum with physically distinct behavior. 
First, at |p| < A, the dispersions of the NG bosons are 
governed by the broken symmetry and match the clas- 
sification into type I and type II. Second, for momenta 
larger than a few times the gap but smaller than the 
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FIG. 7: Coefficients of the type-II NG boson dispersion rela- 
tion at zero temperature, a/m is plotted by the solid line and 
b by the dashed one. 



fcrmion mass m, the dispersions of all NG bosons become 
quadratic. This is in accord with the Bogolyubov theory 
of nonrelativistic BEC [46|. Finally, for |p| > m we enter 
the relativistic regime and the dispersions become linear 
again. (Although it cannot be seen in Fig. [HI the NG 
bosons disappear in the continuum at |p| « 3.6m. So in 
a linear scale, this relativistic domain actually covers a 
large part of the dispersion relations.) 

Let us also remark that in the BCS regime, the fre- 
quency marking the onset of the two-particle continuum 
remains strictly constant for all values of momentum dis- 
played in Fig. [5] The reason is that only in the BCS 
regime has the dispersion of the fermionic quasiparticles 
a nontrivial minimum. Consequently, for all momenta |p| 
up to the Fermi momentum, hp, the pair mode may de- 
cay into two fermions both lying on their Fermi surfaces. 
So the continuum sets at to — A in pn and at lu — 2A in 

P33- 

Another hint on the range of momentum where the 
universal predictions of the broken symmetry apply is 
provided by the size of the coefficients a, b defined by Eq. 
(TTT)) . In Fig. [7] we plot these coefficients against the in- 
verse coupling. Apparently, the coefficient a which gives 
rise to the quadratic onset of the dispersion of the type- 
II NG boson, is suppressed in both the BCS and RBEC 
regimes. In the BCS limit, this is because the popula- 
tions of the gapped and ungapped quarks are almost the 
same and hence the density of the non-Abelian charge 
As goes to zero. In the RBEC limit, the reason is the 
above-mentioned particle-antiparticle symmetry. 



IV. CONCLUSION 

In this paper we studied the relativistic BCS-BEC 
crossover in a class of NJL-type models. Within the 
mean-field approximation, we investigated the superfluid 
phase with particular focus on the collective bosonic exci- 
tations. The formalism we developed is general enough to 
cover all possible spin-zero pairing patterns. The calcu- 
lation of the thermodynamic potential up to one bosonic 
loop splits into a kinematic and an algebraic part. The 
kinematic part is universal and the choice of the pairing 
pattern enters only through the values of the gap param- 
eters in the fermionic quasiparticle dispersion relations. 
The algebraic part bears all necessary information about 
the pairing pattern and can be done at once by means of 
a simple trace. 

Despite the apparent generality of our model La- 
grangian |T]), two extensions are necessary before the 
model can claim phenomenological relevance for relativis- 
tic BCS-BEC crossover. First, it would be desirable 
to have different chemical potentials for the individual 
fermion flavors in order to be able to account for ex- 
ternal (e.g. magnetic) fields, or for the requirement of 
overall charge neutrality. Second, we should also include 
fermion-antifcrmion pairing channels so that we can de- 
scribe the constituent quark masses and the competition 
between the chiral and diquark condensates. 

As a particular application, we described the BCS- 
BEC crossover in the two-flavor color superconductor. 
We investigated in detail the spectrum of NG bosons. 
Their low-momentum behavior was checked to comply 
with the general consequences of spontaneous symmetry 
breaking at finite charge density. Concretely, the five 
broken-symmetry generators correspond to one NG bo- 
son with a linear dispersion relation, and a doublet of 
NG bosons with a quadratic dispersion relation. This 
agrees with the general counting rules for the NG bosons 
in Lorentz-noninvariant systems [iol . l42l . |47| . 

At the very end it should be emphasized once again 
that the mean-field approximation itself is not sufficient 
to describe the crossover accurately. However, the inves- 
tigation of the spectrum of bosonic collective states that 
we performed here, is to be understood as a necessary 
first step towards its systematic improvement. 
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APPENDIX A: SOME USEFUL FORMULAS 
1. Energy projectors 

We use the standard energy projectors A k [3§|, asso- 
ciated to the solutions of the Dirac equation for a free 
fcrmion of mass m, with energy eek, e — ±, 



At-i 



1 H 70(7 • k + m) 



We use the following properties of the projectors, 
# ± /i7o - r7i = 70 ^ ( fc o ± M - e£ k) 



e=± 



7oA k7o = Al k , 7 5 A k75 = A 



Matsubara sums 



— e 
-k" 



1 



(iLO n ) 2 
1 



A 2 (iWr, 



AAB E iJ> 

ei,e 2 



N 



eie 2 



5 2 



f( ei A)f(e 2 B) 



N 



e x A + e 2 B n[e x A + e 2 B) 



where f(x) = \ j{eP x + 1) and n{x) = \/(eP x — 1) are the 
Fermi and Bose distributions, respectively, and f2^r is an 
arbitrary (external) bosonic Matsubara frequency. 



In the main text, we use these two respective formulas 
with a = e£ k , b = A = E^ r , and B — E^ s to define 
the quantities i^/(k, q; ifijy) and J^/(k, q; ifljy). Let us 
note that they possess the following symmetry, 



Computation of the fermion loop contribution to the 
inverse boson propagator, involves the following two 
types of fermionic Matsubara sums (48j , 



TV- 



(iw„) 2 - A 2 (icj n - ifl N ) 2 - B 2 

1 v- (^ + eia)(B + e2fe)/(eiA)/(e 2 B) 
^ iOjv + eiA + e 2 S nfeiA + e 2 B) ' 

ei,e a 



^/(k, q; ifl N ) = j/ r e (q, k; ffljy), 
J^/(k,q;in^) = j/ r e (q,k;iO^). 
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